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Abstract 

We consider wave optics effects in the case of microlensing of an 
extended Gaussian source by a point mass under standard assump- 
tions about the incoherent emission of the source elements. We obtain 
analytical expressions for the power spectrum of the microlensed radi- 
ation, which are convenient to obtain asymptotic formulae in the case 
of a large source. If the source centre, the lens and the observer are 
on a straight line, the power spectrum is presented in a closed form 
in terms of the hypergeometric function. In case of a general location 
of the lens and the source the expression for the power spectrum is 
presented in a form of a series. In the high frequency limit the re- 
sults lead to usual expressions of the geometric optics. Approximate 
expressions for the power spectrum in case of a large source and high 
frequencies are obtained. 



1 Introduction 

Gravitational lens systems (GLS) have an important role for investigations 
of mass distributions in the Universe; in particular, the gravitational mi- 
crolensing is a powerful tool in probing planetary and solar-mass objects 
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in the Milky Way (e.g. [201 ESI 123)- Most of the GLS studies is devoted 
to observable effects that can be described within geometric optics approx- 
imation ([151 IS])- The effects of the physical optics, however, could give 
an independent information about the lensing objects. Also, a detection of 
wavelength-dependent effects in GLS can open the way to entirely new tests 
of the electrodynamics in curved space-time. It has been repeatedly pointed 
out that, though for known GLS the wave effects are small, they may have 
an observational interest in the near future. Different theoretical aspects of 
the physical optics effects in GLS have been studied elsewhere (see a bibli- 
ography in monographs [21 [15] and in papers [91 [UJ [121 UHl El EH [23l [14J . It 
has been repeatedly pointed out that, though for known GLS the wave optics 
effects are small, they may have an observational interest in the near future. 
Recently Heyl [51 [6] turned attention that diffractive effects may be impor- 
tant in case of lensing by planetary masses and planetesimals. A number of 
authors (see, e.g., [HIE] and references therein) pay attention to diffractive 
effects in GLS in connection with the future detection of gravitational waves. 

The wave optics signals of the microlensed radiation from any real source 
are blurred because of its coherence properties; therefore they strongly de- 
pend upon the size of the extended source. The diffractive microlensing of an 
extended source has been studied in the case of caustic crossing events [H [23] , 
in the case of point mass lensing [13]; the mutual coherence of different lensed 
images of an extended source near the caustic has been estimated in a series 
of papers by Mandzhos [HI [T21 [T3] . In order to obtain a power spectrum of 
the microlensed radiation from an extended source, most research involves 
numerical calculations. Nevertheless, it would be desirable to have analytic 
results at least for some simple problems. 

In the present paper we consider a Gaussian source under standard as- 
sumptions about correlation properties of incoherent source elements. We 
obtain a power spectrum of a radiation from the source, which is microlensed 
by one point mass. An observed radiation field is obtained in a standard way 
using the Kirhhoff integral; our new elements are related to the power spec- 
trum. As distinct from earlier results by [H] on this problem, we propose 
an analytical representation of the power spectrum, which is convenient in 
case of a large source radius R as compared to the Einstein radius projected 
onto the source plane (Re, s )- In particular, we propose a closed-form rela- 
tion in terms of the hypergeometric function for the power spectrum when 
the lensing mass is projected onto the centre of the source. Then we use 
this relation in order to study a general case of the source disposition and to 
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obtain approximations for the power spectrum for small Re, s /R- We present 
the first orders of this expansion as well as an expansion in powers of 1/u 2 
in the high-frequency case. 

2 Basic relations 

In this section we formulate general suppositions and relations used below. 
Following [3j [H] and many other papers we use standard considerations of 
the diffraction theory and the gravitational lensing (see, e.g., [21 US]); cor- 
respondingly we neglect polarization effects. Calculations of the field are 
performed in the flat space-time background; this is relevant in case of the 
Galactic systems. However the results can be easily extended to the case of 
extragalactic GLSs after a correct redefinition of cosmological distances. 

Leaving aside the polarization effects we describe the radiation field by 
means of one scalar function (p(t, r). One can assume that in the lens plane 
the light rays gain additional phase shifts corresponding to the gravitational 
time delays [2j [15]; after that the radiation field near the observer can be 
derived using the Kirchhoff 's method. 

Furthermore we work in the Cartesian coordinates; the observer, the lens 
and the source centre are situated in a neighborhood of Z-axis in the planes 
z = 0, z = D d and z = D s respectively. As we neglect polarization, we 
describe the source "current" by a scalar function }{t, r); for our consider- 
ations we can assume that this source lies completely in the source plane: 
r = (y,D s ), y 6 R 2 . We also assume that r) is a stochastic process 
having the correlation properties 

<i(t,r)i(t',r')>=5(y-y')I(t-t',y), (1) 

which take into account that the emission from different points of the source 

is incoherent; < ... > represent an ensemble average. 

~ °° 

For the Fourier transforms r) = -7= / r)e wt dt we have then 

n —00 

<j( W ,r)jV,r') >=5(y-y')5(u-u')I(u,y), (2) 

where I(u,y) = J dt I(t, y)e tuJt is an intensity of an element at point y for 
frequency uj. 

The Helmholtz equation for the Fourier transformed radiation field <p(u, r) 
A(p + u 2 Lp = — Airj (c = 1). 
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We consider solution (p(u), r) of this equation describing the radiation from 
the source plane r' = (r'±,D s ), r'± G R 2 . Let D ds = D s — D d be the 
distance from the lens plane to the source plane (in the flat space). Near 
Z-axis (|r'j_| << Dds) in the lens plane r = (r±,Dd) and under assumption 
that R « D^, where R is a source size, we have 



<p(u},r) 



D, 



ds 



d r'_L j(w,r'_ L )exp 



too 



2D 



ds 



(3) 



In order to calculate the field after passing the lens plane we use the 
well known approximation (e.g., [2]) assuming that the radiation gains an 
additional phase shift uit grav after passing the lens plane r = (Dd, r±). Then 
the transfer function is introduced 



T(u,y) = exp{-iut grav (y)}, 

where t grav (y) is the general relativistic time delay, which is the same for all 
frequencies. The solution just after passing the lens plane is 



D. 



ds 



T(u,y') / d 2 yj(uj,y)exp 



IUO- 



(y' 



2D 



ds 



The radiation comes to an observer at the origin r = 0. The field <p(u), 0) is 
calculated using the Kirhhoff-Sommerfeld method using known field <p s (u, y') 
for r = (y',D d ). 



<p(u,0) = 2 ^. D J rf 2 yVs(w,y')exp 



2Da 



r '2 



(4) 



2m D d 

If the lens be absent, then instead of (p s in this formula we must use ([3 
yielding the same unperturbed solution near the observer. 
The field near the observer is then 



<p(u,0) = 



where 



(w,y) = - / cxi> i 



(x - y) 2 x 2 



>d x . 



(6) 



Using (J2J) we obtain 

< 0)<p*(u/, 0) >= <J(w - u/)P(u;) 
with the power spectrum 

1 ^ 2 



P( W ) 



d yl(w,y) |0(w,y) 



(7) 



3 Central Gaussian source and point lensing 
mass 

Relation (j2J) is written for arbitrary sources and gravitational time delays. 
Further we consider the case of one point microlens at r = (0, Dd) with the 
delay time 

Wy) = 2r 9 HWl/L), r g = 2Gm; (8) 

here m is the microlens mass, L is a dimensional parameter which disappears 
in final calculations; further it is omitted. 

We assume the Gaussian brightness distribution over the source for the 
intensity I from equation (j2J): 



J(w,y,r ) 



f(u) ( (y-r o r 
exp - 



rrR 2 



R 2 



(9) 



here r is the source centre in the source plane; function f(cu) is supposed to 
be the same for all source points; this function determines coherence proper- 
ties of every emitting source element. 
On account of equation (jSJ) we have 



Ti 



exp < iuj 



(x-yj 

2D, 



ds 



2D, 



— 2r g In |x|) 



>gTx 



(10) 



This integral can be calculated exactly [31 121 EH] using the confluent hyper- 
geometric function $ (a, c; x) pQ: 



5 



(u,y) = uj ia e 2D dsT(l - ia) 



l—ia 



$ ia, 1; 



iuD d y 2 * 
2D ds D S: 



(2D ds D d i 



l — ia 



$ ( ia, 1; iay 2 /R 2 E , 



11) 



V D s 

where we used the Kummer transformation [I]; a — cor g , R Es = pr^Z)*] 1 / 2 , 
is the Einstein radius projected onto the source plane, D* = D ds D s /D d ; 

y = |y|- 

For the Gaussian source the power spectrum ([7]) has the form 



P(U,T ) 



d 2 y exp 



(y ~ r o) 
R 2 



,2D d D d J ttR* 
The microlensing effect can be described by the ratio 

? = P(u,t )/P (u), 



Kyj 



(12) 



(13) 



where Po is the power spectrum in the absence of microlensing (a = 0). 

If the source centre is at the origin (ro = 0), the integral ( fl2|) can be 
estimated in terms of hypergeometric functions [15] ; the derivation (that can 
be found in Appendix A) yields T = T (a, a) where 

T (a, a) = exp[2crarctan( / 9)] |T(1 — ia)\ 2 F (ia, —ia; 1 ; (1 + f3 2 )~ 



(14) 

— a/a, a — R 2 Es /R 2 , F(a,b;c;x) is the hypergeometric function. Note 



that yfft = R~ 1 y XD*/n is essentially a ratio of the size of the Fresnel zone 
to the source size. 

For large a » 1 and a ~ 0(1) the argument of the hypergeometric 
function is small and we have To(a, a) ~ 2na. 

For a « 1 that is R E , S « R, and bounded a ~ 0(1) it is convenient to 
use an expansion of the hypergeometric function F(a, b; c; x) near the point 
x = 1 for an integer c. Then ( fl4l) takes on the form 



Tn(a,a) 



exp [2 a arct an (/3)] ^ 1 — a 2 s ^2 

n=0 



T(n + l + ia) 



r(l+«r) 



s w [fc n (cr) - Ins] 
{n + l){n\) 2 



(15) 
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where 

k n (a) = 2^(n + l)-il)(n + l+ia)-il>(n + l-i<T) + —r, (16) 

n + 1 

s = a 2 1 (a 2 + a 2 ) , ip(x) = d\nT(x)/dx. 

We shall write the expansion in a up to terms ~ a 2 and a 2 In a that preserve 
a dependence upon the frequency. In this approximation equation (fT5|) can 
be written as 

T (a,a) = l + 2a + 2a 2 -a 2 [k -2ln(a/a)} . (17) 

The expansion ( f!5|) can be rewritten in the form which is convenient to 
look for asymptotic expansions at large frequencies: 



T (a,er) = exp[2crarctan(a/o-)] (l - V] C n (a) S ^"^1 ) , (18) 

[ ^ (n + l)(n!) 2 J 

where s = a 2 s = a 2 /{I + /3 2 ), 



k n (a) = k n (a) + 2ln(a), C n (a) 



T(n + 1 + ia) 



T(l + ia)a n 



n (i+^h as) 



m=0 



CP 1 



Using asymptotic relations for function ip (e.g., pQ) for large arguments, the 
coefficients k n (a) can be expanded in powers of o~ 2 . Therefore, one can write 



T (a, a) = E a^T^a) + 0(a^ m+1 )), 



n=0 



Up to the terms ~ a 2 we have 



fc n ((r) = fc n (oo) - —j 



n(n + 1) + 
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where k n (oc) = 2ip{n + 1) + (n + 1) 1 ; from ( TE?j) it follows that 

n(n + l)(2n + l) 

C n (cr) = H — h ... 

er^ 



(20) 



(21) 



Then we have the geometric optics limit (er = ur g — > oo) 
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T (a, oo) = T 0) (a) 



Ni f „ ~ a 2n \k„(oc) — 21nal ] 9/v _ , . , . 

exp(2a)] ^ 1 - « 2 E ( n ; i ) (n!)2 J | = 2ae 2 °tf 1 (2a), (22) 



n=0 



where we used a representation for modified Bessel function K\ (see, e.g., [3]; 
formula 8.446). A direct calculation of the amplification coefficient from the 
geometric optics for the point mass lens (Appendix B) is in accordance with 
this expression. 

For the next term of the expansion we have 



a) 



k„(oo) — 2 In a 



n 



^(n + l)(n!) s 



n+l)(2n+l) -a 2 (n + 1) 



+a 2 — n{n+\) — 



4 Non-central source 

Now we proceed to derivation of T in case of a non-central source using 
equations (|T4l)-(|T71). Integration over the angular variable yields 



P(u,t ) = (2D d D ds y 2 2f(oj)R 



>-2 



•exp (-^J J dyyexp 
^ ' o 



V 2 \ T /2yro\ 



(23) 



where r$ = \tq\. Using expansion of modified Bessel function of the 1-st kind 
Jo in the form of Taylor series and the substitution y 1 = tR 2 E s we get 

T(a, cr, r Q ) = a exp + ^o^j l r (l - icr)\ 2 ■ 



oo n / \ 2n °9 



n=0 



71! 



(24) 
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Using differentiation with respect to parameter a and taking into account 
equations ([IT]) . ([II]) this can be represented as 



T(a,a, r ) 



exp 



a n+l /rQ X 2n Qn 



ni 



2 \RJ da n 



T (a,a) 



a: 



(25) 



Obviously this representation is workable when tq/R is not large; this will 
be supposed further. This excludes, e.g., the case of a point source. 

Simple analysis shows that, to obtain the inputs up to the orders of a 2 
and a 2 In a into f[2"5"j) , it is sufficient to use the truncated expression (fT7|) . We 
thus get 

T(a, a, r ) = 1 + 2ae- r o/ fl2 + aV^ 2 . 



where 



r. 



^j[2-* (a)-21n(a/a)]-2^ + 2 l7 (^ 



(26) 



E 

n=2 



X 



n(n — l)n! 



(x - l)[£z(x) - C - lnx] - e x + 1 + 2x, 



Ei(x) is the integral exponent, C is the Euler constant. 

Now we proceed to write an asymptotic expansion for T(a, a, tq) for large 
frequencies. Substitution of ( |20|) into ( 1251) leads to an asymptotic series in 
powers of o~ 2 : 



T(a,a,r ) = £ a^T^a, r ) + 0(a- 2 ( M+1 ); 
with the coefficients 



(27) 



a,r J 



exp 



^ em; 

/ n=0 



71! 



2n 



0" 



T 



( m ). 



a 



a 



(28) 



Some caution is needed when we differentiate term by term asymptotic 
expansion (120]) ; however, this can be easily justified in view of the explicit 
structure of the functions involved. The geometric optics limit yields 



T(«,oo,r ) = T(°)( 



a,r ) 



(29) 
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where T^(a, rp) is given by ( 12 8 p for m = 0. A direct calculation of the 
amplification using standard gravitational lensing theory yields the same 
result (see Appendix B). 

The coefficient T^(a), which describes the first correction to (129]) with 
an accuracy up to ~ a 2 , can be obtained either from (128]) or from (T26|) : 



5 Discussion 

We have obtained analytical expressions for the power spectrum of the ra- 
diation of an extended Gaussian source microlensed by a point mass with 
standard assumptions concerning the incoherent emission of the source el- 
ements. If the source centre, the lensing mass and the observer are on a 
straight line, the power spectrum is given by equation (|14p in terms of the 
hypergeometric function. In the case of the general arrangement the result 
is presented in the form of functional series (j25|) . This representation of the 
power spectrum makes it possible to derive approximations up to any ac- 
curacy in case of a sufficiently small a = (Re, s /R) 2 - The representation is 
effective under condition that the source size is not small as compared to 
the distance to the source centre. The opposite case can be treated by the 
method of [TJ], where an expansion of the magnification around the source 
position is used. Representations (j!4|25p are used to derive lowest orders 
of the asymptotic expressions in cases of the small lens (1261) and the high- 
frequencies (|30|) . We have shown by means of a direct calculation that the 
high frequency limit yields exactly the expressions of the geometric optics. 
As we see from ( TTT]l3"Uj) . the first nontrivial wave optics contribution, which is 
dependent upon the frequency, appears in the terms ~ a 2 (though the first 
geometric optics lensing contribution has the order ~ a). 

At the end of the paper we give some estimates for the radio waveband; 
though much smaller wavelengthes also may be of interest (cf., e.g., [IS])- For 
the effects of the wave optics to be significant, one must have a = r g u ~ 1 
and a ~ 1. For a wavelength ~ 1 — 10 cm the first condition is fulfilled if 
the microlens mass is of the order of 1O _5 M . Such planetary mass objects 
must be common in the Milky Way, and there is a number of observational 
confirmations of exosolar planets, including observations by means of the 




(30) 
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microlensing (see, e.g., |20j). In order that a not be small, the source size R 
must be of the order or less than Re, s - This is true for a wide interval of D* in 
the case of the Milky Way objects; though the probability of microlensing of 
a suitable radio source is small (see, however, j6]). For a distant extragalactic 
source at D s w Dd s ~ 10 3 Mpc microlensed by a planet at Da = 10 kpc, then 
Re, s — lO _2 [M/(lO _5 M )] 1 / 2 pc must be of the order of an inhomogeneity 
scale in the radio waveband. The probability of this microlensing can be 
compared to that of microlensing of quasars by Milky Way stars and stellar 
remnants, and it was pointed out [21] that observation of such events could 
be observationally relevant in the near future. 
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Appendix A. Derivation of equation ( j!4j) in the case of the central 
source. 

From equation (10) we have 



1 2 f2D ds D c 



$ ( ia, 1; W/Re, 



(31) 



For ro = the power spectrum ( 1121) equals to 



PHo) 



2D d D ds J R 2 



]dyy™v(-j?) 



(32) 



where an integration over the angular variable is performed. On account 
of ( !3~Tf3~2l) and using the substitution t = y 2 , we obtain the ratio ( TT3l) corre- 
sponding to r = as follows 
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T (at, a) = e a7T |T(1 - io)\ 2 J dte~ l |$ (ia, 1; iat/a)\ 2 , (33) 

where a = R 2 E i /R 2 . 

The integral in formula fl33|) is a special case of an expression that could 
be obtained by means of formula 6.15.22 of book QlJ). Here, in order to 
estimate (I3"3"j) . we shall perform calculations directly for the integral 

oo 

I(a,X,a',X') — /dte"*$(a,l; At)$(a',l; X't) , (34) 
o 

where < Re(a) < 1, < Re(a') < 1, |A'| + |A| < 1. With this aim we need 
the well-known representations for the hypergeometric function 

1 } 

F (G ' b; 1; x) = r(6)r(i-6) / du (1 " ^ _a ^ 1 ( 1 " ^ b = 
1 00 

=wj^)I dT ll+{1 - z)TrTb ~ 1{1+Tr ' ' (35) 

where we used the Beta-function B(a,b) = T(a)T(b)/T(a + b) and the sub- 
stitution t — t/(1 +t) , and for the confluent hypergeometric function 

1 r 

$(a,l; *) = -— due™u a -\l-u)- a . (36) 

±> a, 1 — a ) J 

We use ([36]) for $ (a, 1; At) and $ (a', 1; A't) in fl34]) and we perform 
integration over t after change of the integration order: 



I (a, A, a', A' 



duu a -\l ~uY a ( dv Va S "' 



B(a, l-a)B(a'l-a') J v 'J 1-Xu-X'v 

Integration over <it>, after substitution i; — >• £: u = r/(l + r) and r = 
£(1 — Xu)/(1 — A — Ait), is easily carried out using an integral representation 
for the Beta-function. 

Then substitution u = r]/(l — X + r/) yields 
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(i - a)~°(i - \')- a ' j J t-'ii+vr'- 1 

l{a,A,a , A ) = — — - / ar] 



B(a,l-a) J li + (i- z ) v ] a 

where z = AA'(1 — A) _1 (l — A') -1 . Therefore, on account of the representation 
035]) we get a final relation 

'^■'^ (i-w-Ir ' (37) 

All calculations are fulfilled in the domain of the parameters a, A, a', A' 
where the integrals involved are convergent. However relation (157)) can be 
analytically continued to a more wide domain which includes the values a = 
ia, a' = —ia, A = ia/a, A' = —ia/a. This yields equation (1141) of the main 
text. 



Appendix B. Gaussian source amplification in geometric optics. 

The amplification of a point source by a point mass lens within the geo- 
metric optics in variables of the source plane is well known [T5J [2] : 

y 2 + 2R% 3 
yjy 2 + AR\ . 



In our case this must be convolved with the Gaussian brightness distri- 
bution Q yielding the total amplification 

A = ^ / d^(y)e-*-^ R2 = f 2 e~% J dyy^y)e^I (^.) 

o 

Taking into account (1551) . using the Taylor expansion of I and substitution 
y 2 = tR\ s we have 

A = A(a,r ) = 
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r 2 OO 
n=0 



-l) n a n+1 fr \ 2n d n 



RJ da r 



dye 



-at 



(39) 



The latter integral can be written [1] in terms of modified Bessel function 
of the second kind K\ yielding 



A(a, 0) = a J dte~ 



-at 



t + 2 

JWTTt 



2ae 2a K 1 (2a) 



(40) 



This is the same as T (a, oo) of fl22|) . 

For an arbitrary source position from (139]) we have 



A(a,r ) 



r 2 OO 
n=0 



-l) n a n+1 fr \ 2n d n 



n! 



RJ da n 



A(a,0) 



a 



(41) 



which is the same as T(a, oo,r ) of (|29|) . 
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